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ON  A  VOLUME  INTEGRAL  EQUATION  USED  IN  SOLVING 
3-D  ELECTROMAGNETIC  INTERIOR  SCATTERING 

PROBLEMS 


Introduction 

Given  an  incident  field  E*"°  generated  from,  for  example,  ein  antenna  or  a 
radeir,  the  problem  is  to  estimate  the  amount  of  electromagnetic  energy  depo¬ 
sition  in  a  nearby  human  object.  While  this  is  a  classiccJ  interior  scattering 
problem,  a  versatile  and  efficient  method  for  solving  it  in  realistic  settings  is 
not  generally  available  [1] .  Some  recent  research  in  this  area  can  be  found  in 
[3,  4,  5]  and  their  references.  Understanding  the  details  of  electromagnetic 
deposition  in  humans  is  essential  for  health  and  safety  exposure  considera¬ 
tions.  Knowledge  of  deposition  in  humans  and  non-human  animals  is  very 
important  to  medical  research  on  the  bioeffects  of  radiation  exposure. 

An  approach  to  solving  interior  scattering  problems  is  the  integral  equa¬ 
tion  approach.  One  of  the  advantages  of  this  approach  is  that  the  boundary 
conditions  of  the  problem  are  automatically  incorporated  into  the  integral 
equation,  a  feature  that  is  especially  useful  in  solving  the  corresponding  in¬ 
verse  problem.  For  simple  problems,  the  equivalence  between  the  classical 
Maxwell’s  equations  approach  and  the  integral  equation  approach  can  easily 
be  demonstrated.  However,  for  general  3-D  problem  where  the  scatterers 
may  have  internal  discontinuity  in  their  material  properties,  this  equivalence 
is  no  longer  obvious.  The  purpose  of  this  paper  is  to  re-examine  the  valid¬ 
ity  and  applicabihty  of  a  formrdation  of  a  volume  integral  equation  approach 
commonly  used  for  solving  interior  scattering  problems.  It  will  be  shown  that 
a  popular  heuristic  derivation  of  a  fundamental  volume  integral  equation  is 
wanting,  as  it  is  mathematically  incomplete.  A  more  complete  analysis  of 
the  volume  integral  equation  formulation  is  given.  As  a  by-product,  the  ap- 
phcabihty  of  this  approach  is  claxified.  It  is  concluded  that  this  approach  can 
be  apphed  to  a  broad  class  of  interior  scattering  problems  including  the  cal¬ 
culation  of  electromagnetic  energy  deposition  in  arbitrary  three-dimensional 
inhomogeneous  objects  commonly  encountered  in  the  biomedical  sciences. 
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Preliminary 

We  assume  a  dielectric  body  of  finite  extent  in  free  space  is  occupying  a  region 
Vb  and  is  radiated  by  an  incident  wave.  The  incident  wave  is  produced  by  a 
source  enclosed  in  another  finite  region  V,  which  we  assume  is  far  away  from 
the  body.  The  free  space  outside  Vb  2ind  will  be  denoted  by  Vo- 

In  each  of  the  three  regions  Vb,  V,,  and  Vo,  the  electric  field, ^  E,  the 
magnetic  field,  H,  the  displacement,  D,  and  the  magnetic  induction  B  satisfy 
the  Maxwell’s  equations.  We  assume  the  fields  are  time  harmonic  with  a  time- 
dependent  factor  where  (3  is  fixed  at  either  1  or  -1.  The  body,  possibly 
nonhomogeneous,  is  assumed  to  be  isotropic,  linear,  and  non-magnetic.  Thus, 
if  E(x,  t)  =  3?{E(x)eJ'^"‘}  and  H(x,  t)  =  Sft{H(x)e^'^‘^*}  then  E  and  H  satisfy: 

V  X  E  -f  jj3u!^o^  =  0 
V  X  H  —  j/3u}eE  =  J 


where 

€  — 

J  =  {Jj,0,J6} 

in  y  =  {^,^0,^6},  respectively.  Here  eo,fJ-o  are  the  free  space  permitivity 
and  permeabihty  respectively.  J  is  the  current  density,  which  when  restricted 
to  V,{Vb)  is  J*(J6)  •  Similarly,  the  permitivity  in  V,{Vb)  is  £,(£&).  Since  we 
assume  the  body  is  hnear,  Jj,  =  crE.  Hence,  with  the  usual  modification  of 
£fc  to  include  the  conductivity,  tr,  we  can  assume,  without  loss  of  generality, 

J6  =  0 

The  incident  wave,  being  the  only  field  present  when  the  body  is  absent, 
satisfies  the  time-harmonic  Maxwell’s  equations: 


V  X  E™ -F  =  0 

V  X  -  j/3a;£"‘’E*"‘’  =  J*”” 


where 

£*’“  =  {es,eo,e^o} 

=  {J.,0,0} 
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in  y  =  {Vg,Vo,Vb},  respectively.  (We  cure  making  the  tacit  assumption  that 
J,  is  not  changed  by  the  presence  of  the  scatterer  (body):  J*”^  =  J  in  y .) 
Defining  the  scattered  field  E*  as  E  —  E*””^  and  similarly  H*  as  H  — 
then  (recall  J^,  =  0) 


V  X  E*  +  =  0  (1) 

V  X  H' -  i/?a;(eE  -  =  0  (2) 

Except  possibly  in  Vg  (which  we  will  not  consider),  this  system  can  be  written 
as 


V  X  E'  +  j/3u;jj.oH‘  =  0  (3) 

V  X  H' -  j/3a;eoE*  =  (4) 

where 

From  the  definitions  of  e  and  e*"°,  it  is  clecur  that  Jg,  vanishes  identically 
everywhere  outside  14,  including  14 . 

Combining  Equations  (3)  and  (4),  we  get 

V  X  V  X  E*  -  fc^E*  =  (5) 


where 

Furthermore,  using  the  vector  identity  V  x  V  x  E*  =  V(V  •  E')  —  V^E*, 
we  also  get  the  equation 

V^E'  +  klE‘  =  R  (6) 

where  A!^(r)  =  e(r)/Xoa;^  and 

R(r)  =  -ik\r)  -  kl)  E(r)  +  V(V  -  E(r)) 

To  solve  any  of  these  equivcilent  equations  for  the  scattered  field  uniquely, 
we  need  to  specify  the  boundary  conditions  it  must  satisfy.  Besides  requiring 
E*  to  satisfy  the  radiation  condition,  we  also  require  the  tangentied  com¬ 
ponents  of  E*  and  H'’  (or  equivalently,  V  x  E*,  since  we  are  assuming  a 
non-magnetic  body)  to  be  continuous  across  Sb,  the  boundciry  of  14: 

(n(r)  X  E'(r))^  =  (n(r)  x  E*(r))_ 

n(r)  X  (V  X  E'(r))^  =  n(r)  x  (V  x  E*(r)) 
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(7) 

(8) 


This  follows  from  the  boundary  conditions  on  the  E  and  H  and  the  fact 
that  E*”^  is  continuous  across  Sb-  lu  the  last  set  of  equations,  quantities 
with  a  positive  (negative)  subscript  denote  values  on  the  exterior  (interior) 
boundary. 


A  Heuristic  Formulation 

In  the  unbounded  free  space,  it  is  well-known  that  the  vector  equation 
V  X  V  X  F  -  klF  =  -j(3wfj,o3  (9) 

has  a  unique  solution  satisfying  the  radiation  condition.  In  fact,  it  is  given 

by 

F(r)  =  (3o  {I  +  ^ VV)  •  g{r,  r')  J(r')  dr'  (10) 

where  V  is  the  support  of  J  and  j3o  =  --j(3ufiQ.  Here  the  function  g  is  the 
free-space  Green  function  for  the  three-dimensional  scalar  wave  equation. 
That  is  to  say  it  satisfies  the  equation 


V^g{r,  r')  +  klg{r,  r')  =  -(5(r  -  r') 

(11) 

The  function  g  is  expHcitly  given  by 

|r-r'| 

H 

^  c 

1 

(12) 

Since  Equation  (5)  has  exactly  the  same  form  as  Equation  (9),  it  is  often 
concluded  therefore  that  the  unique  scattered  field  must  also  be  given  by: 

E*(r)  =  /3c  (i  +  VV)  •  /  ^(r,  r')  Je,(r')  dr'  (13) 

«o 

leading  to  an  integral  equation  for  the  total  field  E: 

E(r)  -  (i  +  ^  V V)  •  ^  g{r,  r')  {k^r')  -  kl)  E(r')  dr'  =  E'”‘=(r)  (14) 

While  this  reasoning  is  plausible,  it  is  mathematically  incomplete,  as  the  issue 
of  boundary  conditions  on  Sb  has  not  been  addressed.  In  the  following,  we 
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will  give  a  fuU  derivation  of  an  integral  equation  for  the  electric  field  in  the 
presence  of  a  finite  number  of  bounded  scatterers  in  free  space.  We  assume 
each  scatterer  is  piecewise  continuously  differentiable  in  e{r).  In  this  deriva¬ 
tion  the  role  of  the  boundeiry  conditions,  both  physical  and  mathematiczil 
(due  to  discontinuity  in  e  within  each  scatterer),  will  be  evident. 


A  Detailed  Formulation 

We  assume  the  entire  3-dimensional  space  can  be  represented  as  the  dis¬ 
joint  union  of  a  finite  number  (^-1-2)  of  sub-regions,  14  (A:  =  0, ... ,  All 

sub-regions  are  bounded  except  one.  14  wiU  be  reserved  for  this  unbounded 
sub-region  and  it  represents  the  free  space  exterior  to  all  the  scatterers  to¬ 
gether  with  one  (assumed  for  simplicity)  sub-region  V^+i  that  contains  the 
source  that  generates  the  incident  field.  Thus, 

k=N+l 

R^=  Vk  (15) 

A:=0 

We  assume  e(r)  is  continuously  differentiable  in  each  sub-region  with  possi¬ 
ble  discontinuity  at  the  boundary.  In  genercd  there  will  be  more  sub-regions 
than  there  are  scatterers,  as  each  scatterer  may  have  to  be  represented  by 
more  than  one  sub-region  on  which  £(r)  is  continuously  differentiable.  This 
formulation  is  quite  general  and  includes  the  important  piecewise  homoge¬ 
neous  case,  which  is  often  used  to  model  human  organs  or  even  the  complete 
human  body. 

We  start  with  a  well-known  vairiant  of  the  Green’s  theorem  involving  a 
vector  field  F(r)  and  a  scalar  field  /(r).  Assuming  both  fields  are  defined  on  a 
connected  region  V  with  a  smooth  boundary  surface  S  and  are  continuously 
differentiable,  then  this  particular  variant  of  the  Green’s  theorem  takes  the 
form: 


J  [/V=*F  -  FVV]  dV  =  [  (V  •  F)/n  -  (n  ■  F) V/ 

-n  X  (V  X  F)/  -  (n  X  F)  X  (V/)  ]  dS  (16) 

Here  n  is  the  unit  vector  normal  to  the  surface  5  of  pointing  into  the 
exterior  of  V.  When  more  specificity  is  required,  this  unit  normal  vector  is 
also  denoted  by  n^. 
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In  each  sub-region  14,  (ife  =  0, . . . ,  N),  the  electric  field  E  satisfies  Equa¬ 
tion  (6)  also.  In  the  source  region  Vn+i,  E  satisfies  Equation  (6)  with  a 
modified  R  (  R-f  )  that  accounts  for  the  source  that  generates  the 

incident  field.  In  any  case,  if  we  let  V  =  V4,  F(r)  =  E(r),  and  /(r)  =  ^(r,r') 
in  Equation  (16),  where  r'  is  any  point  not  on  the  boundary  dVk  of  14,  we 
readily  obtain,  using  Equations  (6)  and  (11),  the  integral  equation 


E(r')  S{r\  14)  = 

-/  R(r)9(r,r')<iK  +  St(r') 

(17) 

where 

%.n  =  {J 

for  r  in  the  interior  of  V 
for  r  in  the  exterior  of  V 

(18) 

and 

S)fe(rO  =  Sd(r';  514, 14)  +  SG(r',  -1;  514, 14)  +  Xfc(r') 

(19) 

and  where 

Sd(r';514,14)  = 

f  (V  ■  E  (r))  g{r,  r')nj*  (r)  dS 

J  aVif  *  * 

SG(r',a(*);514,l4)  =  a(r)  (nj*jr)  •  E^Jr))  Vflr(r,  rV-S' 
Xfc(r')  =  -  /  nJ*(r)x(VxE  (r))^(r,r')d5 

-  W  X  (r))  X  V^(r,r')  dS  (20) 

For  r  on  the  boundary  514,  (r)  =  limr'->r  E(r')  for  r'  in  14-  The  same 

holds  for  a(r).  By  definition,  SG(r',  a;  514, 14)  is  linear  in  a.  Hence, 

£  (*) 

SG(r',  -1;  514, 14)  =  SG(r',  - 1;  514, 14)  +  Y*(r')  (21) 

^Vo 

where 

Y*(r')  =  -SG(r'.  H)  (22) 

^Vo 

Equation  (17)  now  becomes 

E(r')<J(r',14)  =  -  /  R(r)^(r,r')dy 

''Vk 
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+  Sd(r';514,V4)  +  SG(r',^^-l;514,Ffc) 

^Vo 

+  Xfc(r')  +  Y4r')  (23) 

The  following  proposition  allows  one  to  convert  the  two  surface  integrals  Sd 
and  SG  in  Equation  (23)  into  volume  integreds  GdV  and  VGd  defined  as 

GdV(r',a(*);Y)  =  VV  •  /  a(r)E(r)  i?(r,  r') 

JV 

VGd(r';7)  =  ■E(j))g{T,T')dV  (24) 


Proposition  A.  Suppose  a(r)  is  continuously  differentiable  in  a  bounded 
region  V  with  smooth  boundary  dV  such  that  V  •  (o'(r)  +  1)  E(r)  =  0  in  V, 
then 

Sd(r';  dV,  V)  +  SG(r',  a(*);  dV,  V)  =  GdV(r',  a(*);  V)  +  VGd(r';  V)  (25) 


Proof.  This  proposition  basically  follows  from  the  apphcation  of  the 
Divergence  Theorem  and  the  use  of  appropriate  vector  identities  involving 
gradient  and  divergence.  In  particular,  one  can  readily  show 


[  (V-E)hgdS  = 
Jav 

[  a{n-E)VgdS  — 
Jav 


•E)g  +  J^{Vg)V  -EdV 

V'V'  ■  J^aEgdV-  j^(V5) V  •  E  dV 

V'  [  E-VagdV 

Jv 

V  J^(a  +  l)gV  EdV 


(26) 


The  proposition  is  proved  by  showing  the  sum  of  the  last  two  integrals  is  zero 
because  of  the  assumption  on  a. 


Now  summing  Equation  (23)  over  all  sub-regions  T4,  (^  =  0, . . . ,  iV  -|- 1)  and 
using  Proposition  A  and  the  definition  of  R  (appropriately  modified  for  the 
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source  region),  we  readily  obtain  for  r'  not  on  dVk,  {k  =  0, . . . ,  N  +  1) 

EM  =  -  kl)E{r)g{r,T')dV 

VV-/  E(r)g(r,r')JK 

M  «o 

fc=JV+l 

+  (X,(r')  +  Yi(r')) 

A;=0 

+  E"”*(r')  (27) 

where  H  is  the  union  of  all  the  sub-regions  representing  the  scatterers  and 

E^em  -g 

grem(r/)  -  _f  j ,{r))  g{r ,t')  dV 

+  S^+i(r')  -  X^+i(r')  -  Y^+i(r') 

Equation  (27)  is  a  representation  of  any  solution  of  the  Maxwell’s  Equations. 
To  determine  the  unique  solution  to  the  scattering  problem,  we  need  to 
impose  appropriate  boundary  conditions:  Equations  (7)  and  (8)  and  the 
radiation  condition. 


Proposition  B.  If  E  satisfies  the  boundary  conditions  in  Equations  (7) 
and  (8)  and  the  radiation  condition,  then 

A!=JV+1 

E  X/fe(rO  =  0  (28) 

Proof.  The  boundary  dVk  of  each  sub-region  14,  (^  =  0,...,N  -1-1) 

is  made  up  of  n*  >  1  smooth  sub-surfaces  Sk,i  {i  —  1, . . .  ,«*).  Since  each 
sub-surface  Sk,i,  if  it  is  bounded,  is  the  surface  of  or  is  part  of  the  surface  of 
exactly  two  sub-regions,  it  appears  exactly  twice  in  the  sum  in  (28).  Suppose 
a  sub-surface  5  is  a  surface  between  two  sub-regions  14  and  14,  then  its  total 
contribution  Cx  to  the  sum  in  (28)  is 

C7;f  =  -  /  [n^^r)  X  (V  X  E^^  (r))  +  n^=(r)  x  (V  x  E^^  (r))]  p(r,rO 

J  s  ^  ^ 
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However,  since  n^^(r)  =  —  nj*(r),  it  follows  from  the  boundary  conditions 
(7)  and  (8)  that  Cx  vanishes.  For  the  sole  unbounded  surface  (the  bounding 
surface  of  Vo  at  infinity),  its  contribution  to  the  sum  in  (28)  also  vanishes 
due  to  the  radiation  condition. 


Proposition  C.  If  E  satisfies  the  boundary  condition  in  Equation  (8)  then 

*!=JV+1 

E  Yfc(r')  =  0  (29) 

fc=0 


Proof.  As  in  Proposition  B,  suppose  a  sub-surface  5  is  a  surface  between 
two  sub-regions  Vi  and  V2,  then  its  total  contribution  Cy  to  the  sum  in  (29) 
is 


Cy  = 


+ 


-  /[^(fl^'(r).E  (r)) 

Js  evo  ^  ‘ 

^(n^>(r)  ■  (r))]  V3(r,r')  dS 


As  a  consequence  of  Equation  (8)  [2],  we  have 

(’■)  ® 

the  familiar  boundary  condition  on  the  normal  component  of  D.  (We  could 
have  made  use  of  Equation  (30)  directly  without  referencing  its  connection 
with  Equation  (8),  but  this  would  have  obscured  the  fact  that  the  bound¬ 
ary  conditions  on  the  tangenticd  components  of  E*  and  H'  alone  uniquely 
determined  the  field.)  Again,  because  n^*(r)  =  — n^^(r),  Cy  vanishes. 


Because  of  Propositions  B  and  C,  Equation  (27)  becomes 
E(r')  =  f  {k^{r)-kl)E{r)g{T,r')dV 

JVf, 
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E(r)p(r,r')<iy 


+  (31) 

Finally,  when  A:^(r)  =  kl{r),  Equation  (31)  reduces  to  E  =  E''"".  Thus, 
E»-em  _  ginc^  definition.  Hence,  for  r'  not  on  any  boundciry, 

E(r')  =  /  (ft^(r)-fco")E(r)5r(r,r0dF 

+  V'V'-/  E(r)g(r,r')dV 

Jv„  H 

+  E‘”=(r')  (32) 

Which  is  exactly  Equation  (14). 


Conclusion 

We  have  given  in  this  report  a  detculed  derivation  of  a  volume  integral 
equation  that  can  be  used  to  solve  interior  scattering  problems  when  the 
material  property  of  the  scatterer  is  piecewise  continuously  differentiable. 

While  it  is  natural  to  make  the  distinction  between  the  physical  (exter¬ 
nal)  boundary  of  the  scatterer  and  its  mathematical  (internal)  boundaries 
(which  arise  when  there  is  a  discontinuity  of  material  property),  all  physi¬ 
cal  boundaries  are  generally  mathematiced  boundaries  and  vice  versa.  Thus, 
mathematical  boundaries  should  be  treated  with  as  much  care  as  with  phys¬ 
ical  boundaries.  Furthermore,  different  types  of  discontinuities  may  require 
different  treatment.  In  fact,  in  MiiUer’s  classical  text  [2],  he  made  a  careful 
distinction  between  the  case  where  the  material  property  has  a  discontinu¬ 
ity  in  the  derivative  of  e(r)  at  the  boundary  and  the  case  where  it  has  a 
simple  jump  in  e(r)  at  the  boundary.  This  suggests  that  care  should  be  exer¬ 
cised  when  deahng  with  scatterers  with  piecewise  continuously  differentiable 
material  property. 

Indeed  the  volume  integral  equation  we  derived  here  is  based  on  a  careful 
treatment  of  boundary  conditions.  The  derivation  clearly  shows  the  roles 
played  by  boundary  conditions  on  both  the  physical  boundary  and  mathe¬ 
matical  boundary. 
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That  our  formulation  turns  out  to  be  exactly  the  same  as  that  derived 
hemristiccJly  (without  addressing  boundary  conditions)  is  also  reassuring. 
Our  formulation  is  valid  provided  the  observation  point  is  not  on  a  boundary, 
internal  or  otherwise. 

Finally,  with  regard  to  implementation.  Equation  (32)  is  generally  not 
used  directly  because  of  the  difficulty  with  differentiating  inexact  data.  How¬ 
ever,  because  of  the  additivity  of  integration  over  volume,  a  Lippman-Schwinger 
type  equation  can  be  obtained  by  moving  /  -f-  VV  under  the  integral  in 
Equation  (32).  This  will  be  valid  provided  the  observation  point  r  is  not 
on  a  boundary  and  that  proper  account  is  taken  of  the  singularity  of  the 
resulting  Dyadic  Green’s  Function  (in  3-D).  Furthermore,  numericcJ  imple¬ 
mentation  of  this  formulation  will  hkely  be  more  accurate  if  fewer  cells  in  the 
mesh  straddle  surfaces  of  discontinuity.  Hence  a  mesh  of  tetrahedra  would 
be  more  suitable  than  a  mesh  of  rectangular  blocks  if  the  scatterer  is  highly 
inhomogeneous,  as  in  the  case  of  the  human  body. 
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